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Abstract 

Recent work suggests the existence of a non-conventional lowest-lying 
scalar nonet containing the ao(980). Then the ao(1450) and also the Kq(1A30) 
are likely candidates to belong to a conventional p-wave qq nonet. However 
a comparison of their properties with those expected on this basis reveals 
a number of puzzling features. It is pointed out that these puzzles can be 
resolved in a natural and robust way by assuming a "bare" conventional p- 
wave scalar qq nonet to mix with a lighter four quark qqqq scalar nonet to 
form new "physical" states. The essential mechanism is driven by the fact 
that the isospinor is lighter than the isovector in the unmixed qqqq multiplet. 
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I. INTRODUCTION 



The identification and interpretation of tlie low lying scalar mesons are questions of 
great current interest. A variety of approaches and models have been explored M- p2 



In the effective chiral Lagrangian approach from which this paper is motivated, a light 
isoscalar cr (560) in addition to the known light isoscalar /o(980) are needed ^ to produce 
a vrvr scattering amplitude which agrees with experiment. Similarly a light strange k, (900) 
state is needed to understand the experimental ttK amplitude. These three particles were 
postulated to form a nonet, taken together with the known isovector ao(980). Consistency 
of this picture with the properties of the ao(980) as seen in r]' rj-K-K decay and as required 
in vrr/ scattering was checked . The pattern of masses, coupling constants and especially 
the isoscalar mixing angle was observed to be much closer to the one expected from a 
four-quark (qqqq) picture rather than from the conventional two-quark (qq) picture for this 
scalar nonet. The four-quark picture was first proposed by Jaffe in the framework of 
the MIT bag model. Very recent experiments ||2^ on the radiative decays — * Trrj'y and 
7r7r7 have been interpreted p5| , p6| as evidence in favor of the four quark picture of the 
low- lying scalars ao(980) and /o(980). 

Now if one adopts the above picture or, as a matter of fact, any other picture in which an 
unconventional non-gg nonet made of the o"(560), k(900), ao(980) and /o(980) exists, there 
is an interesting puzzle concerning the conventional qq scalar nonet. Such a nonet has an 
interpretation in the constituent quark model as a p-wave excitation and should therefore 
share many characteristics of the other p-wave states (the tensor nonet and two axial vector 
nonets with different charge conjugation properties). To see the puzzling features let us focus 
attention on the experimental scalar candidates with non-trivial isospin quantum numbers 
in the greater than 1 GeV energy range. These are the isovector ao(1450) and the strange 
isospinor ii'Q(1430). According to the Particle Data Group survey |^ (see Table 13.2 on 
page 110), they are the likely candidates for a qq scalar nonet. Then one has the following 
unusual features: 

i) The mass of the ao(1450) (presumably ~ ud) is listed as 1474 ± 19 MeV, about 50 
MeV heavier than the strange K^{14:30) (presumably Kq~^ ~ us) which has a listed mass 
of 1429 ± 6 MeV. Our normal expectation is that the replacement of the d quark in a ud 
composite by an s quark should make the resulting state heavier rather than lighter! 

ii) On comparison with the corresponding members of the p-wave J^^ = 2^+ nonet, 
we see that the qq scalar meson candidates are not lighter; specifically m[ao(1474 ± 19)] > 
m[a2(1318.1 ± 0.7)] and m[K*{U29 ± 6)] ^ m[K*{U32.3 ± 1.3)]. Usually it is expected in 
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the constituent quark model that L • S forces should make the spin particle lighter than 
the corresponding spin 2 particle! This is experimentally evident in the (perhaps too simple) 
cc system where m[xc2{lP)] = 3556.17 ± 0.13 MeV and m[xco(l^)] = 3415.1 ± 0.1 MeV. 

iii) If ao(1450) and Kq{1A30) belong to a conventional nonet their decay widths into 
pseudoscalars should be related. Now, only decay modes into two pseudoscalars have been 
observed for these particles: Kq(1430) nK and ao(1450) Tirj^KK and tt//'. As we will 
see later, SU(3) symmetry predicts 

r [ao (1450)] = 1.5ir [K; (1430)] , (1.1) 

assuming that the total widths are saturated by the decay modes mentioned. On the other 
hand the experimental result is 

r [ao (1450)] = (0.92 ± 0.12) T [K* (1430)] , (1.2) 

which clearly differs from the SU(3) prediction. 

In this note we will show that there exists a model which naturally provides a solution 
to these three problems. This model simply consists of allowing the qq nonet to mix with a 
lighter qqqq nonet. Notice that the isovector in the lighter nonet has a structure udss, with 
two strange quarks. On the other hand the isospinor in the lighter nonet has a structure 
usdd, with only one strange quark. Thus, before mixing the lighter nonet will have the 
isovector heavier than the strange isospinor. This situation is illustrated in Fig. 0, where 
the notation is explained. Details will be given later, but we can easily see how the scheme 
works. The two isovectors mix with each other as do the two isospinors. Since the mixing 
of the two levels repels them this explains point (ii), why the qq scalars appear heavier 
than expected. Similarly the qqqq scalars are pushed down in mass. Point (i), the level 
crossing of the qq isovector and isospinor can be simply understood in the perturbation 
theory approximation: since the ao — ag splitting is smaller than the Kq — Kq splitting the 
"energy denominator" for the isovector mixing will be smaller than the one for the isospinor 
mixing. Hence the isovectors will be more strongly repelled. We must assume that the 
ao — Kq splitting is large enough so that there is no level crossing for the lower mass scalars. 
Finally point (iii), the difference in coupling constants of the ii'g(1430) and the ao(1450), 
can be readily understood from the greater "contamination" of the ao(1450) wave- function 
with the four-quark iso-vector state. 

In our present work we do not discuss the isoscalars of the scalar qq nonet. The reason is 
that the experimental situation is rather fluid at the moment, with many candidates. These 
include the /o(1370) (which may actually correspond to two different states), the /o(1500) 
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FIG. 1. Mixing of two nonets - a', K' , a and K stand respectively for the "phys- 
ical" states ao(1450), -^"0(1430), ao(980) and k(900). Kq and are the unmixed 
isospinor and isovector qqqq states, while K'q and Oq are the corresponding un- 
mixed qq states. 

(which may be a glueball state) and the /j(1710). The present scheme suggests a five-fold 
mixing between the cr(560), the /o(980), two heavier qq isoscalar scalars and a glueball. 



II. MIXING FORMALISM AND MASS SPECTRA 



Our interest is in investigating the mass spectra and later the decay properties of scalar 
mesons which are a mixture of "conventional" qq p-wave states and qqqq states. 

In the quark model the usual qq type scalars are grouped into a nonet, say A^', with 

N'i ~ qaq\ (2.1) 

where a and b are flavor indices and qi,q2,q3 = u,d,s. So, for examples, contains 
one strange quark and one antistrange quark, and contain one strange quark and 
one light antiquark while N'\ and N'2 have no strange quarks. One can also construct 
"multiquark" hadrons, an idea originally discussed by Jaffe In this paper we focus on 
a scalar flavor SU(3) nonet of color SU(3) singlet states with quark composition qqqq. 
Within the context of the MIT bag model, Jaffe showed moreover that the lightest such 
scalar nonet should have a mass less than or in the vicinity of 1 GeV due to the strong 
binding energy of the qqqq configuration arising from hyperfine interactions between the 
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quarks. The 4-quark scalar nonet A^, which transforms in an identical manner to A^' under 
flavor SU(3), can naturally be decomposed (this discussion is a summary of Section II of 
1^) in terms of "dual" flavor quarks (actually diquarks): 

K ~ TaT\ (2.2) 

where 

Ta = eabctq", f'^ = e^'^Qc. (2.3) 

So A^l contains no strange quarks, Nf and contain one strange anti-quark each, while 
Nl, N2 and Nf contain two strange constituents each. As explained in the introduction we 
are not including the experimentally ambiguous isoscalars in our present discussion and so 
the pure qq states in A^' of interest are the isovector and isospinor; their charged components 
are (using the notation of Fig. |T]) : 

ao+ ~ ud, K'q'^ ~ us (2.4) 

and the corresponding members of the qqqq nonet A^ are 

qq^ ~ usds, Kq^ ~ udsd. (2.5) 

By simply considering the strange quark content of these states, and also bearing in mind 
that the qq nonet A^' presumably lies in the same mass range as the p-wave axial and tensor 
meson nonets whereas the bag-model indication is that the qqqq nonet A^ should be less 
than about 1 GeV, we expect an ordering of the masses of these states 

rriKo < ruao < ^a[, < rriK^, (2.6) 

as illustrated in Fig. Suppose initially that the scalar meson nonet A^' is "ideally mixed" 
according to the classic idea of Okubo , applied originally to the lightest vector mesons. 
"Ideal mixing" within the nonet may be defined by the following mass terms of an effective 
Lagrangian density for the qq scalars: 

^'mass = -a'TriN'N') - b'TT{N'N'M). (2.7) 

In fact as discussed in we may define a generalized ideal mixing model for the qqqq nonet 
A^ by the mass terms: 

C^ass = -aTriNN) - bTr{NNM). (2.8) 
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Here Ai is the "spurion matrix" [A4 = diag(l, l,x) where x is the ratio of strange to non- 
strange quark masses in the usual interpretation]. It is worth remarking that although ( p.7| ) 
and ( |2.^ ) are similar in appearance, the difference at the quark level between and N' 
manifests itself through opposite signs of b and b' . This can be seen by noting thatQ 



rriKo 



m. 



ao 



(x 



l)b, 

1)6', 



(2.9) 



where the numerical value of x is around 20.5 p9[. Here b' is taken positive while b is taken 
negative; this agrees with counting the number of constituent strange quarks. 

To see whether a mixing between the nonets A^' and can give states whose properties 
reproduce those of the experimental scalar isovector and isospinor candidates above 1 GeV, 
we consider the simplest invariant term which will induce mixing between and A^', namely 



-7Tr (A^A^') 



(2.10) 



For orientation, we first consider the mixing from the point of view of simple perturbation 
theory applied to two two-state systems. For the isovectors and isospinors we have the 2x2 
mixing matrices 



Ml 



7 



7 m^aL 



and 



Mi 



m 



Ko 7 

7 mj^' 



(2.11) 



At second order in perturbation theory we see that the shifts in the square masses for 
the ao — a'g and Kq — Kq systems have magnitudes 



A. 



7^ 



and A 



K 



mj^i — m 



(2.12) 



respectively. Clearly the ordering of the masses in Eq. ( |2.6| ) implies that A^ > A^' and so 
if the gap between mag and rnxo is sufficiently large relative to the a'^ — K'q mass difference 
we will naturally be able to obtain the level-crossing behavior of Fig. 1. 

Next we proceed to an exact treatment. Invariance of the trace of upon diagonaliza- 
tion implies that 



2,2 2,2 



(2.13) 



* Squared masses are being used since we are working in an effective Lagrangian framework. This 
feature does not play any special role. 
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with an analogous equation holding for the isospinors. Using this condition and the eigen- 
value equation 

« - ^3 (^'^o' - ^9 - 7' = 0, (2.14) 



we solve for the masses of the original unmixed states to obtain: 

1 



2 _ _ 

2 



ml + rua''^ =F V {rua''^ - mlf - 



(2.15) 



Analogous equations follow from the diagonalization of M\. These equations may be read 
as giving for each value 7^ the corresponding masses of the unmixed states which will, upon 
inclusion of C^aassi Isad to the experimentally known physical masses. Reality of the masses 
implies that 47^ < Min{(ma/^ — ■m^)^, [vtik''^ — m'^)'^}. 

We are interested in a scenario where the ordering of the unmixed masses is as in Eq. 
( p.6|) . We find that this can happen provided that m^, — < m^-, — (which holds for 
most of the experimentally allowed range of masses) because in this case the behavior of the 
bare masses is as shown in Fig. 2. Of course the bare and physical masses must coincide 
for 7 = 0. We define jmax'^ = ji'm'i' — "^a)^^ ^^e value of 'jmax'^ depends (since rua is very 
accurately known) on the exact value of m^/. For 7 = 'jmax the 1 = 1 states are maximally 
mixed and the unmixed states are degenerate with square masses equal to ^{ml + ma'"^)- 
We see from Fig. 2 that the choice 7 = 'jmax is expected to result in the largest splitting of 
the bare qq masses rriK^"^ — '^ol^ ■ 

A detailed numerical search shows that the largest value of this splitting is in fact ob- 
tained for 7^ = 7max^ = 0.33 GeV"^ and with the choice of physical masses (within the 
allowed "experimental" range) 

ma = 0.9835 GeV, = 1.455 GeV, 

mx = 0.8750 GeV, mx' = 1-435 GeV. (2.16) 

Fig. ^ shows all the allowed points and their corresponding mass splittings. Notice that mx 
is obtained from the analysis of ttK scattering given in This then yields the following 
masses for the unmixed states: 

rriao = rua'^ = 1.24 GeV, 
mKo = 1-06 GeV , m^^^ = 1.31 GeV. (2.17) 

We see that — m^j^ ^ 65 MeV which is comparable with the analogous splitting of 
the tensor and axial families of order 100 MeV. We also notice that in addition to satisfying 
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2.5 




0.5 I ' ' ' ' ' ' 1 

0.1 0.2 0.3 

y (GeV') 

FIG. 2. Evolution, as a function of 7^, of the bare masses needed to produce the 
physical ones. Of course, the bare and physical masses coincide for 7^ = 0. This 
picture corresponds to the case m^^/^ — m/^^ > rria''^ — nia^, which holds for the 
central "experimental" values. Here the plot is shown for the physical masses 
chosen at the end of Section II. The dot-dashed curve is niKo'^, the solid curve 
is ruao'^, the dotted curve is mx''^ and the dashed curve is nia' ^. 
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FIG. 3. Scan showing the values of the bare mass sphtting m^^^ — rUa'^ resulting 
from different experimentally allowed masses of a, a', K and K' and of 7^. The 
mechanism gives the correct ordering for the approximate range 0.26 < 7^ < 
0.38 GeV^ 

the ordering in Eq. ( |2.6| ) [which can be an explanation for puzzle (i)], we can understand 
puzzle (ii) in this picture since the unmixed qq scalar states are lighter than the analogous 
tensors. Specifically, we have that Tn[aQ] < m[a2(1318.1)] and m[K'Q] < m[i^^^(1432.3 ± 1.3)]. 



III. DECAY WIDTHS 



In this section we address the third puzzle presented in the introduction. The total 
widths of the ao(1450) and the ii"o(1430) are listed in the Review of Particle Physics |]27 



as: 



r*°* [K; (1430)] = 287 ± 23 MeV and F*"* [aq (1450)] = 265 ± 13 MeV 



(3.1) 



which implies the ratio in Eq. (|1.2| ). The only listed decay mode of the -ft'o(1430) is nK 
with a branching fraction of (93 ± 10) % which is close to 100%. On the other hand for the 
ao(1450), the experimental knowledge of the exclusive decay modes is less certain; the nr], 
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KK and ttt]' modes are listed as "seen" without stating any branching fractions. In the 
detailed listings the following ratios are presented: 

ao (1450) ^ iTi^l r[ao(1450) ^Wl , , 

— . ^ -^ = 0.88 ±0.23 and J, " , , ( ^ = 0.35 ±0.16. 3.2) 

r [ao (1450) ^ TXT]] r [ao (1450) ^ vrr^] ^ ' 

In this section we also denote the physical state a' = ao(1450) by and the physical state 
K' = ^^0(1430) by Kq. Despite the uncertainty, for the purpose of our analysis we shall 
assume that the tti], KK and Trrj' modes saturate the ao(1450) decays and that their ratios 
expressed above hold as stated. 

Using isoptopic spin invariance the scalar-pseudoscalar-pseudoscalar trilinear interaction 
terms relevant for these decay channels can be written asQ 

- C = '^"'^ d^,KT ■ a*odf,K + 'ya^nv^l ■ d^izd^r] + Ta^Trr^'ao ■ d^-Kd^-q' 

+ 2^ ^q^Kt ■ d^TzK* + h.c.) . (3.3) 

Hence the perturbative decay width of the ii'Q(1430) is 

r (k; ^ .K) = _ _ ,„,.)^ (3.4) 

where q is the center of mass momentum of the decay products. Analogous expressions 
follow for the ao(1450) partial widths. Thus we have that 

r {a* Tcr]) = 0.00997a5.r,', T (a^ ^ vrr/') = 0.00287aS^v', .3 

r (a* ^ KK) = 0.00707„*x^', T {K* nK) = 0.0U3^k^^^k'. 



Let us initially suppose that the ao(1450) and frQ(1430) are members of a hypothetical 
unmixed scalar qq nonet A^', i.e. 7 = 0. Then their decays into two pseudoscalars are 
presumably described by the interaction 

CM'^r^ = 2A'TT{N'd^cl)d^cl)), (3.6) 

where c/)^ is the matrix of pseudoscalar fields. This pure qq coupling term, when expanded 
into individual fields and compared with (|3.3|) above (where the coupling constant conven- 
tions are defined), leads to the identifications: 



^Derivative coupling is being used because we want our Lagrangian to be a piece of a chiral 
invariant object. See Appendix B of S. 
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Ta^Trr, = -2V2cos9pA' , -fa*^r,' = -2V2sm9pA' , Ja'KK = Ik'-kK = -2A' , (3.7) 



where 9p is the pseudoscalar mixing angle, which we take to be 37° Now if we substitute 
into (|3.5| ) we find the qq SU(3) predictions for the ratios of the total widths: 



r {K*, ^ nK) 
and for the partial ao(1450) widths: 



1.51, (3.8) 



717] 



0.55, ^ = 0.16. (3.9) 

We see that while ( |3.9| ) are just a little below the experimentally allowed ratios ( p. 21 ), 
the ratio ( p.8|) is not consistent with the experimental ratio which follows from (|3.1| ). Thus 
considering the ao(1450) and -ft'o(1430) as members of a pure qq SU(3) nonet does not give 
good agreement with experiment. 

Next we study the predictions for the decay widths of the ao(1450) and the frg(1430) 
in the mixing picture of Section II. In |0 we discussed the general SU(3) fiavor invariant 
coupling of members of a scalar nonet to two pseudoscalars. For the case of the qq scalar 
nonet A^' the most standard form is as taken in (|3.6| ) above. However, for the 4-quark nonet, 
A^, a more natural structure which to a first approximation reproduce s the scalar decay 
pattern is: 

Cnh = Ae'''''edefN%(j)td^4. (3.10) 

For qqqq mesons it seems reasonable that the dominant decays will simply be ones that 
involve a "falling-apart" [^, or rearrangement, of the 4 quarks into two qq mesons. So for 
example since A*"! contains no strange quarks one would expect its decay into KK to be 
suppressed. Indeed the Lagrangian (|3.10|) predicts zero coupling of into KK. 



Upon diagonalization of ( |2.11| ) the physical isospinors are K = k(900) and K' 
^^'0(1430) and we take the mixing convention: 

K (900) \ / cosipK —simjjK 
K* (1430) ) \ sini/jK cosi/jk 

Likewise the isovector mass eigenstates are a = ao(980) and a' = ao(1450) given by 




(3.11) 



ao (980) ] _ / cosV^a -smipa 
flo (1450) J I sini/ja cosifja 




(3.12) 



where the mixing angles are obtained as: 
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tan {2iIj 



K) 



27 



and tan (2ipa) 



27 







m 



(3.13) 



ao 



Now if we take the total trilinear interaction Lagrangian density to be the sum of ( |3.6| ) and 
( p.lOj ) and expand the relevant unmixed isovector and isospinor members of and A^' in 
terms of the physical fields using the mixing convention above, we find that [see ( p.3|) and 
compare with the unmixed case 



'ya'wr) = —2 (^sin'?/'aSin6'pA + V^C0s6'pC0S'?/'aA' 

7a*7rj?' = 2 (^sin'?/'aCos6'pA — A/2sin6'pCOS'?/'a^' 

la'KK = -2 {sinifjaA + COSIpaA') , 

Ik'ttK = -2 (smipxA + cosipxA') . 



(3.14) 



Again we calculate the widths using Eqs. ( ^.5| ) and compare their ratios with experiment. 
It turns out to be helpful to begin by analyzing these ratios in different regions. 
First we consider the limit where |^| is large. In this region, 



r {K* ttK) 
r (aS ^ nr]) 



1.444 



2cos2e, 



cosip 



K 



COSlpa 



(3.15) 



The term in the first bracket is what we obtained above from the couplings in (|3.7| ) and so 
gives the prediction ( |3.8| ) which is smaller than experiment. Recalling ( p.l3| ) we see that 
for mixing angles ^a,'4'K ^ 0, ^ (which is all that is needed in this limit since the relative 

- and so for the opposite sign of 7- may be 



-fO 



sign introduced by considering ipa, "ipK ^ 
absorbed in a change of the relative sign of A and A') we will always have that ipa > 'ipK 
since we are considering the scheme where the ordering of the masses is as in Eq. ( p. 61) . 
Hence the ratio ( |3.15| ) always increases relative to the qq prediction in this limit and so the 
ratio of the total widths ( |3.^ ) will decrease, as required. This behavior is independent of the 
choice of experimental masses as long as they satisfy mK'"^ — rriK'^ > niai"^ — ma^ , which is true 
whenever the mechanism works in order to produce the correct level-crossing behavior for 
the masses. For example, for the illustrative input masses considered at the end of Section 
II we have that ipa = \ and iI)k ~ 31° which implies that 

r (a* ^ vrr/) 



T{K*,^7iK) 



0.606 



(3.16) 



giving that 



1.036. 



(3.17) 
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FIG. 4. Plot of ^ against ^. The dot-dashed curve is V [fsT* (1430) -> TrfsT], the 
sohd curve is F [oq (1450) vrr/], the dashed curve is F (1450) — > KK and 
the dotted curve is F [oq (1450) vrr^']. 



Within the errors quoted in this agrees with the experimental result ( p.l|) and is much 
closer than the qq prediction of Eq. ( |3.(j| ). Furthermore from ( pj.l4D for large |^| we have 
the same prediction as in Eq. ( pl9|) . 



Next we suppose conversely that |^| is large. In this region, 



F {Kl ^ -kK) 
F (a5 ^ 



1.444 



cos^6'„ 



sin?/) 



K 



sinT/iQ 



2.115, 



(3.18) 



for il}^ = 31°. The ratios of the ao(1450) widths now become 



T{a* 



KK) 



F {a*o ^ TTT]) 



0.7071- 



sin 9r. 



1.95, 



Ffa^ 



TTT] 



F (aS ^ vrr/) 



0.2828cot2^p = 0.49. (3.19) 



In this limit, where it is the qqqq decay modes of the ao(1450) and the ^^0(1430) that 



dominate, we see that in particular the first ratio in ( |3.19| ) is well outside the experimentally 
allowed range. 

For A' ~ A a graphical analysis is helpful since in this region the ratios of the widths 
blow up. In Fig. 3 we plot the widths themselves (up to an overall normalization of ^). It 
is seen that, for the central values of ( ^TTD and (|3.2|) the correct width order F {Kq nK) > 



Via: 



— Trr^) > F (c 



KK] >T (a* 



TTT] 



is obtained for the "asymptotic" regions ^ > 



1.2 and ^ < 



-3.2. Inside, where —3.2 < ^ < 1-2, the correct width order cannot be 



obtained for the central values. 
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In summary, the above analysis shows that for large |4-| the mechanism significantly 



improves [see for example ( |3.17| )] the ratio / , compared with the prediction based on 
a pure qq description of the ao(1450) and the ii'Q(1430). Outside of this asymptotic region 
more detailed analysis is needed and requires additional experimental guidance. 



IV. NUMERICAL ANALYSIS OF THE DECAY WIDTHS 

In this section we give a more detailed numerical analysis of the decay widths. We will 
take into account the experimental uncertainties for comparison with theory. Furthermore 
we will include a more general form of decay interaction. Finally the decay widths of the 
lighter scalars K = fi;(900) and a = ao(980) will also be discussed. The input masses will be 
kept the same as in Eq. ( |2.16D , and we will continue to use 7^ = 7max^ = 0.33 GeV^. The 



general interaction Lagrangian describing the decay widths has the form (see [Q) 

£,„i. = Ae'''''edefN%(f)ld^<pf + CTT{Nd^(f))TT{d,(f)) 

+ A'e''''edefN'%<Pld^<Pi + C"Tr(iV'9^0)Tr(9^0) + ■ • • (4.1) 

where the three dots stand for terms which do not contribute to isovector or isospinor decays. 

We first consider the limit C = and C = 2A', as in Section III. In this limit, the above 
Lagrangian simplifies | to 

Lint. = Ae'''"'edefNtd^ct>ld^(t>i + 2A'Tr(iV'a^09^0) + ■ • • (4.2) 

We scan the AA' parameter space numerically and search for regions consistent with the 
available experimental data on the decay widths of these scalars. We start with the ao(1450) 
decay widths as they impose the strongest restrictions on the parameter space. First, we 
find that the experimental estimate r*°*' [ao(1450)] = 265 ± 13 MeV restricts A and A' to the 
perimeter of the ellipse shown in Fig. ^ - the thickness of the perimeter is related to the 13 
MeV uncertainty of the decay width. We then search for regions that are consistent with 
the current experimental estimates on the ratios r[ao(1450) — > -ft'-ft']/r[ao(1450) Tirj] = 
0.88±0.23, and r[ao(1450) -> 7rV]/r[ao(1450) ttt]] = 0.35±0.16. Regions consistent with 
the first and second ratios are respectively shown by dark and light shading. The vertical 



•tjt is helpful to use the identity: 

A'e''''''edefN'%r,d^4 = 2A'Tv{N'd^^d^cl)) - A'TT{N')TT{d^^d^cl)) - 2^'Tr(iV'5^</.)Tr(9^,/.) + 
A'Tr(iV')T¥(5^0)Tr(9^</>). 
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axis corresponds to the conventional interaction term for qq nonets, whereas the horizontal 
axis represents a natural interaction for qqqq nonets, as previously discussed. We see in 
Fig. ^ that within our model we cannot exactly describe the current experimental data on 
the partial decay widths of ao(1450). Obviously a natural four-quark interaction is far from 
the allowed regions, while a natural two-quark interaction seems to be a favorable scenario 
for description of the available experimental data. Also small distortions from the natural 
two-quark interaction, although slightly improving the situation, do not exactly describe 
the data. This is described more quantitatively in Table |. We have fitted the prediction of 
our model for the total decay width as well as the decay ratios, to the above experimental 
estimates and searched for the best values of A and A' . The natural two-quark interaction 
(column one) is compared with the more general case that natural four-quark interactions 
are also allowed (column two). Although the of the fit gets slightly reduced, effectively 
the best point remains around the natural two-quark interaction. 

In order to see whether we could get a better description of ao(1450), we have also 
extended our investigations to the more general case where C 7^ and C 7^ 2^4', i.e. working 
with the general Lagrangian ( |4.1| ). The result is given in column three of Table |, indicating 
that even with the introduction of more general interaction terms, the current experimental 
data is still not exactly described. We also notice that in this case C — 2A' ^ and that 
\2A'\ dominates A and C. Thus the best fit in this case also is similar to column I, although 
the fit is slightly improved. Therefore, the simplified model Eq. ( |4.2| ) already provides 
a reasonable picture for understanding the nature of the ao(1450). We should emphasize, 
however, that the available experimental estimates of the decay channels of ao(1450) are not 
very accurately known. More accurate experimental data on ao(1450) would be useful for 
our purposes. 

Next, we include the i^^^(1430) in the picture. We take experimental values F*"* = 
287 ± 10 ± 21, and F^i^/F*°* = 93 ± 10%, and search for regions that give T[K^{U30) 
ttK] ^ 267 ±50 MeV. These are shown in Fig. ^ with two parallel strips in the north-west to 
south-east direction. We see in the figure that within our model ( |4.2| ) there are overlaps of 
regions in parameter space AA' that explain most of the decay properties of both ao(1450) 
and i^'o*(1430). 

Now that we can understand the decay widths of the heavier scalars, we explore the 
possibility of explaining the decay widths of the light scalars within the same theoretical 
setup. We proceed by further exploring the parameter space AA' in the limit C = and 
C = 2A', for regions that explain decay properties of the lighter physical nonet members 
ao(980) and /t(900). We search for regions consistent with F[ao(980) 711]] ^ Q5±5 MeV in 
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A(GeV ') 

FIG. 5. Regions in the A A' parameter space [see Eq. ( [4 .21) ] consistent with 
the currently available experimental estimates on the decay widths of ao(1450). 
Points on the ellipse are consistent with the total decay width of ao(1450). 
Dark and light gray regions respectively represent points consistent with the 
experimental ratio r[ao(1450) KK]/T[ao{U50) ttt]] = 0.88 ± 0.23, and 
r[ao(1450) ^ 7rr/']/r[ao(1450) ^ irr]] = 0.35 ± 0.16. 
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Fitted Parameters 


Eq. 


Eq. i3) 


Eq. (4.1) 







0.10 ±0.12 


1.03 ±0.12 




-2.55 ±0.06 


-2.60 ±0.06 


-3.53 ±0.12 










1.36 ±0.27 




2A' 


2A' 


-6.56 ±0.27 


Predicted Decay Widths 


r*°*[ao(1450)](MeV) 


265 


265 


265 


r[ao(1450) ^ KK]/T[ao{U50) nr]] 


0.55 


0.53 


0.53 


r[ao(1450) ^ 7rr/']/r[ao(1450) nr]] 


0.16 


0.18 


0.18 




1.161 


1.157 


1.157 



TABLE I. Best numerical values for the free parameters in the scalar pseu- 
doscalar pseudoscalar interaction Lagrangian, found by fitting the prediction of 
our model for the total decay width and ratio of the partial decay widths of 
ao(1450) to the experimental data. The first column corresponds to an inter- 
action natural for qq, while in the second column interaction terms natural for 
qqqq are also included. In the third column the more general interaction Eq. 
(4.1) is considered. 
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agreement with experimental measurement in pO| as well as the theoretical estimate in [Q. 
We also search for regions consistent with ^ r[fi;(900) — t^K] ^ 40 ± 5 MeV in agreement 
with theoretical estimates of the properties of «:(900) given in |^. The result is also shown in 
Fig. 1^, indicating that there are regions in the parameter space of our model (A ^ ±1 and 
A' q=3) that are approximately consistent with the decay properties of the light scalars in 
addition to those of the heavy scalars. 

We have given in Table ^ our best fits for A, A', C and C, resulting from comparing 
our theoretical prediction to the experimental data. We have also displayed in the same 
table the predicted decay widths. In the limit C = and C = 2 A' (column one), the 
resulting decay widths have the right order of magnitude, although some of them are not 
within the ranges allowed by experiment. The fit gives ^ = —2.4 and so, as expected from 
the discussion of Fig. ^ the widths F (1450) — > KK and F [ao (1450) nrj] have the 
wrong order. Outside this limit (column two), we get a better agreement with experiment 
(as of fit ^-Iso indicates), and except for the ratio F[ao(1450) i^'-ft']/F[ao(1450) nr]], 
all other decay widths are within their experimentally allowed ranges. We notice that in the 
general case (column two), C ~ 2A', which means that the decay interaction for nonet A^' 
remains close to that natural for qq. We also notice that in this general case, C ^ 0, which 
is expected from our previous results on decays of the low lying light scalars 0. 

V. DISCUSSION 

(i) We studied the properties of the ao(1450) and i^'o(1430) scalar mesons (which are 
usually considered to belong to a conventional p-wave qq nonet in the quark model) in a 
framework where a lighter scalar nonet (of qqqq type) was also present. It was found that 
certain puzzling features of these two particles could be naturally explained if the qq and 
qqqq nonets mix with each other to form new physical states. The essential mechanism is 
driven simply by the fact that the isospinor is lighter than the isovector in the unmixed qqqq 
multiplet. 

(ii) Although we carried out the analysis in a qqqq picture for the unmixed light scalar 
nonet, it seems reasonable that it could also be done for other models of the light scalars (like 
the unitarized quark model P,P|, or molecular models |^) in which they have somewhat 



^This is a width corresponding to the numerator, rather than denominator of a partial wave 
ampUtude as explained in Q. 
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-5-3-1 1 3 5 

A(GeV ') 

FIG. 6. Regions in the A A' parameter space [see Eq. (|4.2|) ] consistent with the 
current experimental and theoretical estimates on the decay widths of ao(1450), 
i^^^(1430), ao(980) and k(900). Points on the ellipse are consistent with the 
total decay width of ao(1450). Squares and circles respectively represent points 
consistent with the experimental ratio r[ao(1450) KK]/r[aQ{14:50) ttt]] = 
0.88 ± 0.23, and r[ao(1450) ^ 7rr7']/r[ao(1450) ^ vrr/] = 0.35 ± 0.16. 
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Fitted Parameters 


Eq. (13) 


Eq. (0) 




1.40 ±0.12 


1.19 ±0.16 




-3.26 ±0.07 


-3.37 ±0.16 


C(GeV-^) 





1.05 ±0.49 


C"(GeV-^) 


2A' 


-6.87 ±0.50 


Predicted Decay Widths 


r*°*[ao(1450)](MeV) 


274 


263 


r[ao(1450) ^ KK]/T[aoiU50) ^ vrr/] 


0.30 


0.42 


r[ao(1450) ^ 7rr/']/r[ao(1450) ^ vrr^] 


0.52 


0.32 


T[K*{U30) nK]{MeV) 


245 


298 


r[ao(980) ^ nri]{MeV) 


57 


65 


r[/t(900) ^ 7rir](MeV) 


45 


41 


r*"* [ao( 1450)] /r[ir*( 1430) ttK] 
(expected value : 0.99 ± 0.24) 


1.12 


0.88 




1.864 


0.757 



TABLE II. Best numerical values for the free parameters in the scalar- 
pseudoscalar-pseudoscalar interaction Lagrangian, found by fitting the predic- 
tion of our model for both the low lying and next-highest scalars to the experi- 
mental data. The first and second columns correspond respectively to the limit 
(C = and C = 2 A') and to the general case outside this limit. 
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different four-quarlc interpretations. 

(iii) We did not investigate the heavier isoscalar particles because the experimental sit- 
uation is still rather ambiguous. Clearly this is an interesting future project. 

(iv) In our treatment we used the simplest mixing term ( |2.10 ) and obtained fairly good 



agreement with experiment. The model can easily be generalized to include different mixing 
terms in the effective Lagrangian; for example 

Tr [A^ (A^A^' + A^'A^)] , Tr (A^) Tr (A^') , (5-1) 
Tr (MN) Tr (A^') , Tr (A^A^') Tr (A^) . 

(v) Although our focus in this paper was on the heavier scalars, the model of course 
describes the lighter ones too. If we want to describe only the lighter scalars, as in , we 
can imagine "integrating out" the heavier scalars. In the simplest approximation, based on 
neglecting the symmetry breaking terms in Eqs. (p.7|) and (p.8|), we would just replace 



A^' — > -—N. (5.2) 

To check the consistency of this with our previous work we might ask how much the decay 
coefficient of the light scalars [A in Eq. ( |3.1CI| )] gets modified due to this replacement. Using 



]6D and the identity in footnote 3 then gives 



A — ^ A - —A'. (5.3) 
2a' ^ ' 

Using A = 1.2 GeV"-^ and A' = -3.4 GeV~^ from column 2 of Table II, together with 
72 = 0.33 GeV^ and a' = 0.76 GeV^ from dgj) shows that A = 1.2 GeV"^ in the present 
paper is to be replaced by A = 2.5 GeV~^ in a model where the heavy scalars have been 
eliminated. This is in rough agreement with A = 2.9 GeV~^ found in Table III of 
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